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Fertilization of sea urchin eggs involves an increase in protein synthesis associated with a 
decrease in the amount of the translation initiation inhibitor 4E-BR A highly simple reaction 
model for the regulation of protein synthesis was built and was used to simulate the 
physiological changes in the total 4E-BP amount observed during time after fertilization. 
Our study evidenced that two changes occurring at fertilization are necessary to fit with 
experimental data. The first change was an 8-fold increase in the dissociation parameter 
(koffi) of the elF4E:4E-BP complex. The second was an important 32.5-fold activation 
of the degradation mechanism of the protein 4E-BR Additionally, the changes in both 
processes should occur in 5min time interval post-fertilization. To validate the model, 
we checked that the kinetic of the predicted 4.2-fold increase of elF4E:elF4G complex 
concentration at fertilization matched the increase of protein synthesis experimentally 
observed after fertilization (6.6-fold, SD = 2.3, n = 8). The minimal model was also 
used to simulate changes observed after fertilization in the presence of rapamycin, a 
FRAP/mTOR inhibitor. The model showed that the elF4E:4E-BP complex destabilization 
was impacted and surprisingly, that the mechanism of 4E-BP degradation was also 
strongly affected, therefore suggesting that both processes are controlled by the protein 
kinase FRAR/mTOR. 

Keywords: translational control, sea urchin embryos, mechanisms of fertilization, deterministic model, translation 
simulation 



INTRODUCTION 

Early development of the sea urchin embryo is one among the 
models which contribute to establish the paradigms of the mech- 
anisms of translation, the eukaryotic universal process for protein 
biosynthesis (Mathews et al., 2007). Many factors are involved in 
all steps of translation machinery, namely initiation, elongation, 
and termination and have been extensively reviewed (Mathews 
et al., 2007). For the great majority of the mRNAs, the initia- 
tion step of translation occurs on activated mRNAs that contain 
a m''GpppN molecule (where N is any of the four nucleotides) 
at the 5' end, also known as a cap structure, thus referred 
as cap-dependent translation initiation (Jackson et al., 1995). 
Not exhaustively, protein translation initiation begins with the 
binding of the eukaryotic initiation factor 4E (eIF4E) to the 



cap-structure of the mRNA. The protein eIF4E recruits a large 
scaffolding protein called eIF4G that interacts, among others, 
with eIF4A and eIF3 Hnking the 5' end of the mRNA and the 43S 
preinitiation complex (Gingras et al., 1999a). After scanning from 
the 5' end to the start codon, the first amino acid is incorporated 
and the peptide elongation step proceeds (Jackson et al., 2010, 
2012; Hinnebusch, 2011). The protein eIF4E is a major target 
for regulation of translation initiation (Sonenberg and Gingras, 
1998). Its availability in the cells depends on the presence of 
4E-BPs, eIF4E binding proteins, which compete with eIF4G for a 
common binding site on eIF4E (Mader et al., 1995). Thus 4E-BPs 
sequester eIF4E, and consequently prevent cap-dependent trans- 
lation (Haghighat et al., 1995). The 4E-BPs phosphorylation sta- 
tus regulates their interaction with eIF4E: underphosphorylated 
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4E-BPs bind to eIF4E and inhibit cap-dependent translation, 
whereas hyperphosphorylated forms do not (Pause et al., 1994). 
Most of the data available to date indicate that FRAP/mTOR 
(FKBP12 and Rapamycin-Associated Protein/mammalian Target 
Of Rapamycin) is the main kinase that phosphorylates 4E-BP 
(Brunn et al., 1997; Burnett et al., 1998) on four conserved 
residues in metazoans (Gingras et al., 1999b). 

Among the experimental advantages of the sea urchin embryo 
model, genes encoding all factors involved in translation are 
present and are non-redundant; thus, while mammalian genomes 
contain three 4E-BP homologs, sea urchin contains an unique 
form (Morales et al, 2006). At fertilization, an increase in protein 
synthesis occurs, which is necessary for the initiation of devel- 
opment and for the occurrence of the first cell cycle (Epel, 1990; 
Cormier et al, 2003; Gilbert, 2003). The protein synthesis increase 
is independent of new transcription and involves cap-dependent 
initiation from maternal mRNAs already present in the unfer- 
tilized egg (Cormier et al., 2003). Fertilization triggers 4E-BP 
release from eIF4E and consequently protein synthesis transla- 
tion through a rapamycin-, and thus, mTOR-sensitive pathway 
(Cormier et al, 2001; Salaun et al., 2003, 2004). An original 
mechanism of 4E-BP regulation first demonstrated in sea urchin, 
corresponds to the rapid disappearance of 4E-BP pool following 
fertilization (Salaun et al, 2003, 2004, 2005; Oulhen et al, 2007, 
2009). Although fertilization is known to provoke many changes 
in the egg metabolism that could interfere on translation reg- 
ulation (Epel, 1990; Gilbert, 2003; Parrington et al, 2007), we 
demonstrate here that a reduced model involving 4E-BP, eIF4E, 
and eIF4G is sufficient to explain the changes observed at fertiliza- 
tion on 4E-BP level and on the increase of cap-dependent trans- 
lation. Simulations using this minimal model suggest that both 
a 32.5-fold increase of the 4E-BP degradation mechanism and 
an 8-fold change in the dissociation constant (koffi ) of eIF4E:4E- 
BP complex are required at fertilization and should occur in a 
5 min time interval after fertilization. Furthermore, simulations 
of changes after fertilization in the presence of rapamycin show 
that the drug strongly affects the eIF4E:4E-BP complex stability 
and, unexpectedly, the 4E-BP degradation mechanism implying 
that both mechanisms are controlled by the FRAP/mTOR protein 
kinase. 

MATERIALS AND METHODS 
PROTEIN PRODUCTION AND PURIFICATION 

The construction of the plasmid pAr encoding eIF4El from Mus 
musculus has already been described (Pyronnet et al., 1999) and 
the protocol to produce and purify Flag-tagged MmIF4El was 
performed as described herein. The cloning of Strongylocentrotus 
purpuratus 4E-BP and the production of the recombinant His- 
tagged Sp4E-BP protein was performed as described in Gosselin 
et al. (2011). Briefly, recombinant Sp4E-BP was produced in BL21 
(DE3) strain E.coli (Novagen) and purified on an affinity column, 
packed with a chelating sepharose Fast Flow Resin preloaded with 
Ni^+ ions (Amersham Pharmacia Biotech). A second purifica- 
tion step was performed on a Superdex 75 column (Amersham 
Pharmacia Biotech). After concentration with Amicon® Ultra-4 
Centrifugal Filter Units (lOkDa), concentration of the purified 
proteins was estimated with a Nanodrop ND-1000. 



Purification of native eIF4E was performed from 
Sphaerechinus granulans egg extracts using m''GTP beads as 
described previously (Oulhen et al., 2010). After washing, beads 
were suspended in Laemmli loading buffer. The amount of 
recovered S^IF4E was determined after resolution on SDS- 
PAGE electrophoresis and revelation by Coomassie blue staining. 
Quantification was performed by comparison with a range of BSA 
concentrations migrated on the same gel. The bands were quan- 
tified using the ImageJ 1.43] program (Wayne Rasband, National 
Institutes of Health, USA) after digitization of the stained gel. 

PREPARATION OF GAMETES AND FERTILIZATION 

Sphaerechinus granularis sea urchins collected in the Brest area 
(France) were maintained in running seawater. Spawning of 
gametes, fertilization, and cell culture were as described (Oulhen 
et al., 2010). When indicated, 20 |iM rapamycin (LC laborato- 
ries) from a 20 mM stock solution in ethanol was added to the 
eggs 20 min before fertilization. Cleavage was scored under a light 
microscope. Each experiment used gametes from a single female 
exhibiting greater than 90% fertilization. 

QUANTIFICATION OF 4E-BP AND elF4E IN SEA URCHIN EGGS 

To quantify eIF4E and 4E-BP in the sea urchin unfertilized eggs, 
total crude protein extracts were prepared by direct dissolution 
of 20 [i\ egg pellet (corresponding to 12,000 eggs) with 150 [xl of 
SDS-Fix buffer as described (Belle et al., 201 1). A range of differ- 
ent extract volumes was submitted to SDS-PAGE and Western- 
blotted together with known amounts of purified SgIF4E or 
recombinant Sp4E-BP as prepared above. Proteins in the different 
samples were immuno- revealed using mouse monoclonal anti- 
body directed against eIF4E from rabbit (BD transduction labo- 
ratories, Lexington, KY) or rabbit polyclonal antibodies directed 
against human 4E-BP2 (Rousseau et al., 1996), a generous gift 
from Nahum Sonenberg (McGill University, Montreal, Quebec, 
Canada). The antigen-antibody complex was measured by chemi- 
luminescence using horseradish peroxydase-coupled secondary 
antibodies according to the manufacturer's instructions (ECL or 
ECL-F, GE Healthcare Life Sciences). Quantification of the bands 
was done using the ImageJ 1.43j program. Measurement has been 
performed on extracts from eggs of 9 different females for eIF4E 
and 4 females for 4E-BP. Protein concentrations were calculated 
considering an egg volume of 0.3 nl (diameter 80 |im). 

DETERMINATION OF PROTEIN SYNTHESIS IN VIVO 

Unfertilized eggs (5% suspension in sea water) were incubated 
in the presence of [^^S]- L- methionine during 1 h at the final 
concentration of 10(jLCi/ml. When indicated, 20|xM rapamycin 
was added to the incubation mixture 20 min prior to fertiliza- 
tion. Eggs were then harvested by centrifugation, rinsed three 
times, re-suspended in Millipore-filtered sea water in the pres- 
ence or absence of rapamycin and fertilized. At different times, 
batches (500 of the embryo suspension were taken; the cells 
were pelleted and frozen in liquid nitrogen. Soluble protein 
extracts were prepared as described in Oulhen et al. (2010). 
[^^S]- L-methionine incorporation into proteins was measured 
on duplicate aliquots of the extracts after 10% TCA precipitation 
on Whatman 3 M filters and counting in a scintillation counter in 
the presence of Optiphase Supermbc scintillation liquid. 
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MEASUREMENT OF 4E-BP LEVEL AFTER FERTILIZATION 

At different times, following fertilization, 20 [il pelleted embryos 
were dissolved in 150 |il of SDS-Fix buffer to prepare total crude 
protein extracts as above. Proteins were separated on SDS-PAGE 
and 4E-BP was analyzed by immunoblotting using the rabbit 
polyclonal antibodies directed against human 4E-BP (see above). 
The antigen-antibody complex was revealed by chemolumines- 
cence, the signal was digitalized and quantified using the ImageJ 
1.43j software. Densitometric analysis gave relative numerical val- 
ues, which were used to calculate percentages of 4E-BP present in 
embryos considering as a reference (i.e., 100%) the quantity of 
total 4E-BP in unfertilized eggs. When necessary, percents were 
converted into absolute concentrations based on values for 4E-BP 
calculated as described above. 

SURFACE PLASMON RESONANCE (SPR) DATA ACOUISITION 

Data were acquired with Reichert SR7000DC spectrometer 
instrument (Buffalo, USA). Running buffer was HEPES 50 mM 
pH 7.7 KCl 150 mM 0.5% Igepal and flow rate was 25 (jiL.min^^ 
60 (jlRIU of Sp4E-BP was immobilized on a mixed SAM (1 Cn- 
(OEG)6-COOH:10 Cii-(OEG)3-OH) (Buffalo, USA) via clas- 
sical amine coupling chemistry. For interaction experiment, 
concentration of Flag-MmIF4El:m''GTP was assayed by spec- 
trophotometry and stock solution was centrifuged before use. 
Then a concentration range from 4.37 |xM to 53.9 nM of Flag- 
MmIF4El:m''GTP (3-fold dilution series) was injected on the 
Sp4E-BP chip at 25°C. Following each Flag-MOTlF4El:m^GTP 
injection, the chip was regenerated with 30 s of 10 mM NaOH. 
Each curve displayed was double referenced with a set of 6 
blank buffer injections. Data were then fitted using a global 
analysis method with Scrubber 2.0a software (Biologic Software, 
Australia). Fitting errors are reported for each kinetic constant. 

RESULTS 

THE MINIMAL MODEL ELABORATION 

A model for cap-dependent translation initiation in sea urchin 
has been proposed (Belle et al, 2010) based on the systemic 
environment of BIOCHAM software (Calzone et al, 2006). The 
Biocham model contains too many reactions and undetermined 
parameters to allow quantitative simulations of dynamic changes 
occurring at fertilization. We therefore investigated the strategy 
of extracting a minimal core model. The model was constructed 
with the following reactions: (Rl) the reversible association 
between eIF4E and 4E-BP, (R2) the reversible association of eIF4E 
to eIF4G, (R3) the irreversible production of protein by the ini- 
tiation complex containing associated eIF4E:eIF4G, (R4) the syn- 
thesis of 4E-BP by the initiation complex containing associated 
eIF4E:eIF4G and (R5) the degradation of free 4E-BP 

(Rl) eIF4E -I- 4E-BP <-> eIF4E:4E-BP (koni/koffi = 1/KDi) 
(R2) eIF4E -|- eIF4G <-> eIF4E:eIF4G (kon2/koff2 = I/KD2) 
(R3) eIF4E:eIF4G -> Protein -|- eIF4E:eIF4G (kcatProtein) 
(R4) eIF4E:eIF4G -> 4E-BP -|- eIF4E:eIF4G (kcat4£-Bp) 
(R5) 4E-BP -> 0 (kiys4E-Bp) 

In this model, eIF4E is considered to be always associated to 
mRNAs since its affinity for the cap structure and the stability 



of the eIF4E: cap mRNA complex were demonstrated to be very 
high in other species (Sonenberg and Gingras, 1998). Global pro- 
tein production is assumed to be proportional to the amount 
of eIF4E:eIF4G complex (R4). The total amounts of eIF4E and 
eIF4G were shown to remain constant after fertilization in sea 
urchin (Oulhen et al, 2007). 

The differential equations resulting from the model using the 
law of mass action are as follow: 

(1) d[eIF4E]/dt = kom [eIF4E:4E-BP] -|- koffz [eIF4E:eIF4G] - 
koni [eIF4E] [4E-BP] - kon2[eIF4E] [eIF4G] 

(2) d[4E-BP]/dt = koffi [eIF4E:4E-BP] -|- kcat4E-BP [eIF4E:eIF4G] 
- ktysis4E-Bp[4E-BP] - koni [eIF4E] [4E-BP] 

(3) d[eIF4G]/dt = koff2[eIF4E:eIF4G] - kon2[eIF4E][eIF4G] 

(4) d[eIF4E:4E-BP]/dt = ko„i[eIF4E][4E-BP] - k^ffi [eIF4E:4E- 
BP] 

(5) d[eIF4E:eIF4G]/dt = kon2[eIF4E] [eIF4G] 
koff2[eIF4E:eIF4G] 

(6) d[Protein]/dt = kcatProtein[eIF4E:eIF4G] 

PARAMETER DETERMINATION IN UNFERTILIZED EGGS 

The relevance of the core model in unfertilized eggs was stud- 
ied by combining experimental data acquisition and learning 
procedure of unknown parameters. Initial total concentrations 
of 4E-BP, eIF4E, and eIF4G as well as kinetic parameters of 
the association of eIF4E:4E-BP and eIF4G:eIF4E complexes had 
to be experimentally determined and were further used to fit 
parameters for 4E-BP synthesis and degradation. 

Concentrations of the translation factors 

The concentrations of 4E-BP and eIF4E in unfertilized eggs were 
determined by immunoblotting and densitometric quantification 
as indicated in the Materials and Methods section. The total con- 
centration of 4E-BP and eIF4E were respectively 3.67 \iM (SD = 
0.83, n = 4), and 2.15 |xM (SD = 0.29, « = 6) thus indicating 
an excess of total 4E-BP over total eIF4E. The concentration of 
eIF4G could not be determined using the procedure described 
for 4E-BP and eIF4E for two reasons: first, sea urchin eIF4G 
was present under complex post-translational modified forms 
(Oulhen et al., 2007) and it was therefore not clear which bands 
could be ascribed to eIF4G, and second, no significant and mea- 
surable amounts of endogenous nor recombinant protein had 
been obtained. However, it was previously demonstrated that 
after fertilization a rapid exchange of 4E-BP for eIF4G occurred 
on eIF4E with a stoichiometric behavior (Oulhen et al, 2007). 
Therefore, simulations were performed using a minimal eIF4G 
concentration equivalent to eIF4E concentration 2.15 [lM. 

Determination oftlie Idnetic parameters 

The kinetic parameters for 4E-BP: eIF4E interaction were deter- 
mined by surface plasmon resonance (SPR). 4E-BP was immo- 
bilized on a chip and several concentrations of eIF4E were 
injected onto the chip. Figure 1 shows the overlapped sensor- 
grams in relative Index of Refraction Unit (RIU) and their 
associated theoretical fits (solid lines) resulting from global 
data analysis according to a one step 1:1 reaction. The kinetic 
parameters determined were koni = 9.3 ± 0.1 x 10^ M^^s^\ 
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koffi = 2.2 ±0.1 xlO^"* s^' and the affinity constant KDi = 
23 nM (Table 1), all in good concordance with published data in 
other species (Abiko et al., 2007; Mizuno et al., 2008). Regarding 
the high conservation of the binding sites among species (Joshi 
et al, 2005) the measured parameters likely reflect the values in 
the sea urchin Sphaerechinus granulans. 

A direct measurement of eIF4E:eIF4G interaction could not 
be analyzed by surface plasmon resonance. A kon2 = 1-82 x 10^ 



M^i.s"' and a koff2 = 2.0 x 10"'' s"' for eIF4E:eIF4G interac- 
tion were used (Table 1) from published data in other species 
(Von Der Haar et al., 2006; Umenaga et al., 2011), resulting in 
determining an affinity constant for KD2 = 110 nM (Table 1). 

Parameters for 4E-BP synthesis and degradation 

The constants for the synthesis of 4E-BP (kcat4E-Bp) and its 
degradation (kiys4E_Bp) had to be determined. The degradation 
constant of reaction R5 was calculated from experiments consist- 
ing in the measurement of the total amount of 4E-BP remaining 
after treatment of unfertilized eggs in the presence of the protein 
synthesis inhibitor emetine. In that condition, the disappearance 
of 4E-BP depends on reactions Rl, R2, and R5 since no pro- 
duction could occur (R4). The disappearance of total 4E-BP was 
simulated and the best fitting value for kiys4£-BP was determined 
by optimizing the least square error to the experimental values 
(Figure 2). The best fit was obtained for a value of 5.9 x 10"* 
s '(Table 1). 

Since 4E-BP remains constant in the untreated unfertilized 
egg, the synthesis constant, kcat4E-Bp(R4), could be calculated 
using reactions Rl to R5. A least-square method was used to con- 
straint a constant value, as close as possible to 100%, for the total 
amount of 4E-BP over time; the inferred catalytic constant for 
4E-BP synthesis was kcat4E-BP = 3.2 x 10"^ s"' (Table 1). 

Calculation of factor concentrations 

The unfertilized eggs are in a steady state in which the total 
concentrations of 4E-BP, eIF4E, and eIF4G remain constant as 
verified from immunoblotting experiments over the period of 
conservation of the eggs. Using reactions Rl and R2, the concen- 
trations of free and complexed forms of 4E-BP, eIF4E, and eIF4G 
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FIGURE 1 I Surface plasmon resonance Interaction of elF4E and 4E-BR 

Each sensorgram, relative refractive index (m-RIU) as a function of time, was 
obtained after injection of elF4E from 4.37 ixM to 53.9 nM according to a 
three-fold dilution range. Overlapped signals were fitted using a single 
equilibrium kinetic isotherm (solid red lines). 



Table 1 | Parameters of the model for translation changes at fertilization (see text for details). 



Unfertilized eggs Fertilized Fertilized eggs Method of determination 

eggs + rapamycin 





9.3 X 10^ M-'s"' 






SPR 


koffi 


2.2 X 10"^ s"'' 


x8 


x1 


SPR and simulation 


KDi 


23 nM 








kon2 


1.82 X 103 M-is-i 






SPR-Bibliography 


koff2 


2.0 X IQ-^s"'' 






SPR-Bibliography 


KD2 


IIOnM 








klys4E-BP 


5.9 X 10-" s-'' 


x32.5 


x16 


Emetine treatment of eggs and simulation 


kcat4E-BP 


3.2 X 10-3 s-i 






Simulation 


Time interval of 




5 min 


33 min 


Simulation 


kinetic parameter 










change 










4E-BPto,3i 


3.67 ixM 






Immunoblots densitometry 


elF4E total 


2.15iiM 






Immunoblots densitometry 


elF4Gtotai 


2.15iiM 






Bibliography/Estimation 


4E-BP 


1.90 M.M 


x0.11 




Simulation/model perturbation 


elF4E 


0.02 iiM 


x12.9 




Simulation/model perturbation 


elF4G 


1.79 liM 


xO.36 




Simulation/model perturbation 


elF4E:4E-BP 


1.77 liM 


x0.20 




Simulation/model perturbation 


elF4E:elF4G 


0.36 M-M 


x4.2 


x2.4 


Simulation/model perturbation 


Protein synthesis 


1 


x6.6 


x3.5 


Radioactive rate of methionine incorporation 
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FIGURE 2 I Determination of thie degradation rate of 4E-BP in 
unfertilized eggs. Unfertilized eggs were incubated in the presence of 
100 |iM emetine. The amount of total 4E-BP remaining upon time was 
determined from immunoblotting experiments as indicated in the material 
and methods section. The values from 8 independent experiments are 
plotted as black cross. Simulation using the model (reactions Rl, R2, and 
R5) leading to the best fit is shown as red line. The total amount of 4E-BP is 
expressed as % of initial value. 



could be calculated. From the differential reactions (4) and (5) of 
the model: 

(4) d[eIF4E:4E-BP]/dt = koni [eIF4E] [4E-BP] - koffi [eIF4E:4E- 
BP] =0 

(5) d[eIF4E:eIF4G]dt = kon2[eIF4E][eIF4G] - koff2[eIF4E: 
eIF4G] = 0, 

and knowing the total concentrations: 

[4E-BP],ot,i = 3.67 (iM, [eIF4E],ot,i = [eIF4G],otal = 2.15 [iM, 
the two equations become: 

koni ( [eIF4E] total- [eIF4E:4E-BP] - [eIF4E:eIF4G] ) ( [4E-BP] 
total- [eIF4E:4E-BP])-koffi[eIF4E:4E-BP] = 0 
kon2 ( [eIF4E] total- [eIF4E:4E-BP] - [eIF4E:eIF4G] ) ( [eIF4G] 
total- [eIF4E:eIF4G])-koff2[eIF4E:eIF4G] = 0 

from which the unknown values [eIF4E:4E-BP] and 
[eIF4E:eIF4G] can be calculated. From these, the free amounts of 
[eIF4E], [4E-BP],and [eIF4G]are deduced from the total amount 
minus the complexed forms 

[eIF4E] = [eIF4E]total - [eIF4E:4E-BP]- [eIF4E:eIF4G] 
[4E-BP] = [4E-BP] total - [eIF4E:4E-BP] 
[eIF4G] = [eIF4G] total - [eIF4E:eIF4G] 

The calculated values are all reported in Table 1. 

At that point, the parameters and steady state concentrations 
of all the reactions of the model are known in the unfertil- 
ized eggs (Table 1). Since protein synthesis was determined by 
[^^S]-methionine incorporation into all translated proteins and 



expressed as a percentage of incorporation, the kcatProtein could 
not be determined. However, since in the model, the amount 
of protein synthesized is proportional to the concentration of 
[eIF4E:eIF4G], changes in [eIF4E:eIF4G] will reflect the changes 
in the protein synthesis rates. 

LEARNING PARAMETERS AND CONCENTRATIONS AFTER 
FERTILIZATION 

The kinetic of 4E-BP disappearance observed in vivo reveals two 
main features: a decrease of the total 4E-BP protein starting from 
fertilization and lasting about 15 min, followed by a stabiliza- 
tion of the level of the protein at 18-20% of the initial amount 
(Figure 3A). In parallel, the rate of protein synthesis increases and 
stabilizes around 30 min after fertilization resulting in the linear 
accumulation of neosynthetized proteins (Figure 3B). Therefore, 
fertilization corresponds to a time -dependent change from a 
steady state (unfertilized state) to another one (fertilized state) 
reached after 30 min. With respect to the minimal model pro- 
posed, two non-exclusive changes could explain the experimental 
values of total 4E-BP evolution: a destabilization of eIF4E:4E-BP 
complex (Rl) and an increase in the degradation rate kiysis4E^BP 
of 4E-BP. The two hypotheses were simulated separately. Since 
the biochemical modifications induced by fertilization are not 
instantaneous (Epel, 1990; Cormier et al., 2003; Gilbert, 2003), 
the time necessary for the biochemical modifications to occur at 
fertilization through the transduction pathway was introduced in 
the model as a parameter time change varying linearly from the 
initial value of parameters in unfertilized eggs to their final value 
after fertilization. 

A set of model simulations was first performed by per- 
turbing the value koffi (Rl) ranging from its unfertilized egg 
value (Table 1) up to 100 times the value and by introduc- 
ing a parameter time change ranging between 1 and 15 min. 
For each simulation, the mean-square error between the model 
prediction and 4E-BP experimental data was computed. The 
least- worst fit (sum of square residual = 46712) is shown in 
Figure 4A and, not surprisingly, the simulation led to a non- 
biologically relevant increase in the total amount of 4E-BP. 
A second set of simulations was performed by perturbing the 
value of kiysis4E-Bpranging from its value in unfertilized eggs 
(Table 1) to 100 times the value while still introducing a param- 
eter time change ranging from 1 and 15 min. The least- worst fit 
(sum of square residual = 3122) is shown in Figure 4B and the 
corresponding simulation did not fully match the experimental 
data. 

Therefore, simulations were performed by perturbing both the 
koffi value (Rl), and the kiysis4E-BP (R4)values together with a 
parameter time change ranging from 1 to 15 min. Using least 
square method, a very good fit with the experimental data was 
obtained (sum of square residual = 1207). The best fit was 
obtained when koffi increased 8-fold from the value in unfertil- 
ized eggs and kiysis4£-BP increased 32.5-fold with a parameter time 
change of 5 min (Figure 4C). Figure 4D shows in a color scale the 
least square sum of square residual values obtained in the simu- 
lations for the ranging values of koffi and kiysis4E-Bp. The best fit 
(white cross) was located in a small area related to the accuracy of 
the determinations. 
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FIGURE 3 I Changes in 4E-BP total amount and protein synthesis in vivo 
after fertilization. (A) 4E-BP was measured by immunoblotting and 
densitometric quantification as described in the material and methods 
section. Mean values from 10 independent experiments are plotted as a 
function of time. The total amount of 4E-BP is expressed as % of initial value. 



(B) Protein synthesis was determined from incorporation into proteins of 
I'^^Slmethionine at different times after fertilization as indicated in the 
Materials and Methods section. Mean values of protein accumulation from 11 
independent determinations are plotted as a function of time in arbitrary 
unit (UA). 
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FIGURE 4 I Simulation of fertilization changes on 4E-BP using the 
minimal model. Best fits (red curves) compared to experimental data 
(black cross) are shown: (A) Ranging koffi from unfertilized value (1 to 
100-fold) and the parameter time change (1 to 15 min). (B) Ranging 
k|ys4E-BP from unfertilized value (1-100-fold) and the parameter time 
change (1-15 min). (C) Changing koffi and k|y54E-BP from unfertilized value 
(1-100-fold) and the parameter time change (l-15min). (D) Color 



representation of the fit according to the change of k^ffi and k|ys4E-Bp 
when the parameter time change is fixed to 5 min (kcat4E-BP and koffi 
ratio are the ratio between the value reached after fertilization and the 
unfertilized value). The color scale on the right was associated to 
distance to data (sum of square residual) obtained in the simulations. 
The white cross corresponds to the best fit. The total amount of 4E-BP 
is expressed as % of initial value. 



An 8-fold increase in koffi leads to a 8-fold increase in the 
KD1/KD2 ratio (Rl and R2), such an increase would also be 
obtained through a kon2 8-fold increase or a koni or koffi 8-fold 
decrease. We therefore performed simulations with these different 
parameters, associated with the 32.5-fold increase in kiysis4E-Bp 
and a parameter time change of 5 min. The best fits (sum of 
square residual) for each parameter were respectively 1207 for 
koffi, 3692 for kon2> 3419 for koni and 4214 for koff2 (Figure 5) 



indicating that only the 8-fold increase in koffi led to a good sum 
of square residual value. 

Figure 6 shows the distance to data (sum of square residual 
value and SD) calculated when each parameter varied by steps 
from 1/32 to 32-fold compared to unfertilized values. Simulation 
with the best fit (sum of square residual = 1189) was obtained 
with respectively a koffi x 8, koni x 1, and kon2 and koff2 x 1/27. 
Regarding that this sum of square residual is not significantly 
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FIGURE 5 I Simulations performed varying the indicated parameters 

(kpffi, kon2. l<oni. l<off2) 3* fertilization. k|ys4E-Bp was fixed as 32.5-fold 
increase compared to unfertilized eggs. Parameter time change from 



unfertilized values was set at 5 min. The best fits (red line) for each parameter 
are shown (A-D) as compared to experimental data (black cross). The total 
amount of 4E-BP is expressed as % of initial value. 



different from the distance (1207) obtained with the unique 
koffi variation, the koffi role in eIF4E:4E-BP destabUization at 
fertilization appears to be the most relevant. 

Altogether, the model allows to conclude that destabilization 
of eIF4E:4E-BP complex by an 8-fold increase of its dissoci- 
ation rate (koffi) associated with a 32.5-fold increase in the 
mechanism of 4E-BP degradation, both taking place in 5 min fol- 
lowing fertilization, are necessary and sufficient to explain the 
observed experimental changes occurring in the total amount 
of4E-BP. 

Using the minimal model, the parameters determined in the 
unfertilized eggs and the changes predicted in the fertilized eggs 
(Table 1) the kinetic of concentration changes of all parame- 
ters were simulated (Figure 7). Importantly, the concentration of 
eIF4E:eIF4G complex increased 4.2-fold from its initial value in 
the unfertilized eggs, in a time interval of 30 min. The result- 
ing calculated increase and the protein accumulation kinetic 
(Figure 7, black line) were highly compatible with the in vivo 
protein synthesis changes after fertilization (see Figure 3B). The 
4.2-fold increase in eIF4E:eIF4G complex compares favorably 
with the experimental 6.6-fold increase in protein synthesis mea- 
sure as described in the material and methods section (SD = 
2.3, n = 8 independent experiments). Although we are aware that 
reaction R4 of the model is a very highly simplified representa- 
tion of the mechanism of cap-dependent translation, and that 
many other changes on translation factors have been reported 
at fertilization in sea urchin (Oulhen et al, 2007; Belle et al., 
2011; Costache et al., 2012) the assumption that translation was 



proportional to the amount of eIF4E:eIF4G complex appeared 
rational in the minimal model elaboration. 

SIMULATION OF RAPAMYCIN EFFECT ON 4E-BP CHANGES AND 
PROTEIN SYNTHESIS INCREASE AT FERTILIZATION 

It has been already reported that FRAP/mTOR is involved in the 
protein synthesis increase following fertilization (Salaun et al, 
2003, 2004, 2005; Oulhen et al, 2007). The effect of rapamycin, an 
inhibitor of the protein kinase FRAP/mTOR, has been analyzed 
in early development of sea urchin embryo. Protein synthesis 
measured as described in the Materials and Methods section was 
53.0% (SD = 5.1, M = 8) of the value in control fertilized eggs. 
Rapamycin slowed down the decrease in the total amount 4E-BP, 
which stabilized at 41.7% (SD = 5.4, n = 8) of unfertilized value 
after 30 min. This observation was rather surprising since this 
value is higher than the value in fertilized eggs (18%), while pro- 
tein synthesis was decreased as compared to control fertilized 
eggs. 

The changes observed in the presence of rapamycin were mod- 
eled as a modification of the parameters of the minimal model. 
Simulations were performed with koffi ranging from 1 to 8-fold, 
kiysis4E-BP from 1 to 32.5-fold and changes for time parameter 
from 1 to 100 min. A good fit was obtained for a kofti 1-fold the 
value in unfertilized eggs, kiys;s4E-BP value 16-fold the unfertil- 
ized value and a parameter time change of 33 min (Figure 8 and 
Table 1). In these conditions, the concentration of eIF4E:eIF4G 
at steady state was 2.58-fold compared to the 4.25-fold in fer- 
tilized eggs in the absence of rapamycin. The model therefore 
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FIGURE 6 I Distance to data calculations from simulations performed 
varying together all of the parameters k^ff^, k^nz, k„„i, k^m- Each 
indicated parameter varied by steps from 1/32 to 32-fold compared to 
unfertilized values in such a way that a 8-fold increase in KD1/KD2 was 
maintained constant (this gives 10434 possibilities). The k|ys4E-Bp was 



fixed as 32.5-fold increase compared to unfertilized eggs and the 
parameter time change was fixed to 5min. For each possible ratio 
(fertilized value vs. unfertilized value), the mean and SD of the 
distances to data were reported. Green bar corresponds to region 
where simulation fits well. 



predicts a rapamycin-induced limited increase in protein synthe- 
sis of 53.6%, in very good concordance with the experimental data 
value of 52.0%. 

Therefore the model predicted that not only koffi increase 
after fertilization but also the degradation mechanism of 4E-BP 
is under the control of FRAP/mTOR pathway. 

DISCUSSION 

While the regulation of protein synthesis is a highly complex 
process, we show here that using a minimal model based on 
three main actors eIF4E, eIF4G, and 4E-BP, the changes occur- 
ring at fertilization in sea urchin could be fully simulated with 
adjustment of a minimum of parameters. The model postulates 
that fertilization triggers two events: (1) an increase in the KDi 
of the eIF4E:4E-BP complex due to an 8-fold increase of the 
dissociation constant (koffi) and (2) a 32.5-fold increase in the 
mechanism of 4E-BP degradation (kiysis4£-Bp). Therefore, the dis- 
appearance of 4E-BP cannot result from the unique dissociation 
of the eIF4E:4E-BP complex. Interestingly, in mammalian cells it 
has been found that 4E-BP1 is degraded in eIF4E-knock-down 
cells (Yanagiya et al, 2012). The authors show that the non- 
eIF4E-bound hyperphosphorylated 4E-BP1 is degraded while 
eIF4E-bound hypophosphorylated 4E-BP1 is stable reinforcing 
our assumption in the model that only free 4E-BP is degraded. 
However, it should be noted that there is no indication in this 
study of an activation of the mechanism of degradation itself 
while it is the case in sea urchin fertilization. 



As expected, the model predicts that rapamycin reverses the 
fertilization-induced change in kgffi. An attractive hypothesis 
would be that fertilization in sea urchin activates the kinase 
FRAP/mTOR that phosphorylates 4E-BP and increases the koffi 
of its association with eIF4E, resulting in a destabilization of 
eIF4E:4E-BP complex. A role for FRAP/mTOR protein kinase on 
the destabilization of eIF4E:4E-BP complex at fertilization is com- 
patible with the observations ( 1 ) that rapamycin, a FRAP/mTOR 
inhibitor, inhibits protein synthesis increase in sea urchin fertil- 
ization, and inhibits the disappearance of 4E-BP upon fertiliza- 
tion and (2) that the 4E-BP protein was shown to be a substrate 
for FRAP/mTOR in other organisms; in mammals, multiple 
and hierarchical phosphorylation events provoke the release of 
4E-BP1 from eIF4E, in concordance with a KDi increase (Gingras 
et al, 2001). Although the relationship between FRAP/mTOR 
activity and the affinity of 4E-BP for eIF4E is not established 
in sea urchin embryos (Oulhen et al., 2009), it has been shown 
that the degradation of 4E-BP correlates with its phosphorylation 
(Cormier et al., 2001). Therefore, the possibility that fertiliza- 
tion activates FRAP/mTOR resulting in a koffi increase should 
be experimentally analyzed when the purified proteins will be 
available under their non-phosphorylated and phosphorylated 
forms. 

We cannot exclude that fertilization could have also mod- 
ified KD2. Since eIF4G undergoes electrophoretic changes 
at fertilization (Oulhen et al, 2007) it would be of inter- 
est to identify the post-translational changes occurring 
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FIGURE 7 I Simulation of the concentration changes occurring at 
fertilization. The concentration changes of the indicated constituents were 
calculated and plotted vs. time after fertilization using the parameters 
depicted in Table 1; "elF4E," "elF4G," and "4E-BP" are the concentrations 
of the free forms of these proteins; "elF4E:elF4G" and "elF4E:4E-BP" 
correspond to the concentrations of the complexed forms; "Total 4E-BP" 
represents the total concentration of 4E-BP (free and complexed). 
Normalization was done with 100% corresponding to the initial 
concentration of total 4E-BP i.e., 3.67 \^M (see Table 1). The kinetic of 
accumulation of protein was calculated from the kinetic of elF4E:4E-BP 
changes and is shown in relative units (black curve). 
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FIGURE 8 I Simulation of fertilization changes on 4E-BP in the 
presence of rapamycin. Simulated curve (red line) obtained with optimized 
parameters (see text) compared to experimental data (black cross). The 
total amount of 4E-BP is expressed as % of initial value. 



on this factor, and its consequences on eIF4E:eIF4G 
interaction. 

Furthermore and unexpectedly, the model also predicts that 
the mechanism of 4E-BP degradation was also rapamycin- 
sensitive. Very interestingly, the degradation mechanism of 4E-BP 



appears to be a FRAP/mTOR target, which could also be bio- 
chemically tested when the mechanism will be identified. For 
now, this feature is described for the first time using sea urchin 
fertilization and would be of great interest to be searched in other 
species and other physiological regulations. 
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